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In this paper we study breakage rate statistics of small colloidal aggregates in non-homogeneous 
anisotropic turbulence. We use pseudo-spectral direct numerical simulation of turbulent channel 
flow and Lagrangian tracking to follow the motion of the aggregates, modelled as sub-Kolmogorov 
massless particles. We focus specifically on the effects produced by ductile rupture: This rupture is 
initially activated when fluctuating hydrodynamic stresses exceed a critical value, a > a^r, and is 
brought to completion when the energy absorbed by the aggregate meets the critical breakage value. 
We show that ductile rupture breakage rates are significantly reduced with respect to the case of 
instantaneous brittle rupture (i.e. breakage occurs as soon as ct > (Jcr). These discrepancies are due 
to the different energy values at play as well as to the statistical features of energy distribution in 
the anisotropic turbulence case examined. 


I. INTRODUCTION 

Breakage rates of micro and nano aggregates in turbulent flow are of high relevance to a wide variety of applications. 
These range from industrial processes, such as operations involving colloids [1], compounding of plastic and elastomeric 
materials [2], and dispersion of ceramics [3], to environmental processes, e.g. sedimentation of marine snow [4] or 
formation of marine aggregates [5]. In almost all cases of practical interest, breakage is caused by two mechanisms. 
The first is usually referred to as impact breakage and is caused by energetic collisions of the aggregates with other 
aggregates or walls. This mechanism, however, is not relevant to the present study because it applies to aggregates 
that are large with respect to the characteristic length scale of the fluid shear (the Kolmogorov length scale if the flow 
is turbulent) and have a significant density difference compared to the fluid. The second mechanism is usually referred 
to as shear breakage and applies to small aggregates with density close to that of the fluid. In this case, breakage is 
caused by aggregate deformation induced by the hydrodynamic stresses. Such deformation generates internal stresses 
that can break the aggregate following two processes: If the response time of the aggregate to deformation is very 
small then breakage is instantaneous, and the aggregate can be referred to as brittle; if breakage depends on the stress 
history and a significant amount of energy is required to overcome deformation, then the aggregate can be referred 
to as ductile. 

In either brittle or ductile rupture, the phenomenology of turbulent breakage is still not fully understood because the 
complexity of the flow field adds to the intricacy of the aggregate morphology in determining how the hydrodynamic 
forces redistribute over the structure of the aggregate and how stresses accumulate in critical locations where the 
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cohesive force that keeps the primary particles of the aggregate together can be overcome. To provide a basic 
understanding of turbulent breakage, many investigations (see for instance the recent works by [6-9] and the references 
therein) have focused on the influence that the hydrodynamic stresses have on the rate at which breakage occurs, 
neglecting the details of the aggregate inner structure. This choice [10] was motivated by the observation that 
the occurrence, in the surroundings of a small aggregate, of instantaneous stresses capable of inducing breakage 
is controlled by the statistical (spatial and temporal) distribution of these stresses, which is in turn dictated by 
the turbulent flow field [7]. However, in the size range of interest for the present study (aggregates smaller than 
the Kolmogorov length), analyses were carried out considering brittle aggregates and instantaneous breakage [6, 8]. 
This assumption would be fully justified in highly viscous flows, such as dispersions in liquid polymers, where the 
stresses required to break the aggregate are very low [11]. But in low-viscosity systems the effective hydrodynamic 
stress required for breakage changes significantly depending on the nature of the flow and rupture can frequently 
be determined by the stress history [12], namely by what we refer to as ductile rupture effects in this paper. As a 
result, current knowledge of the breakage process for the case of sub-Kolmogorov ductile aggregates has remained 
hitherto elusive in spite of its practical importance in areas such as industrial materials processing [13] and rheology 
of dense colloids [14]. Design of such processes usually relies on semi-empirical correlations based on idealized flow 
configurations [15-19]. When turbulence is present, it is generally modelled with a single value of the kinetic energy 
dissipation, a rather crude assumption already for homogeneous and isotropic flows. However, real facilities involve 
mixing in containers, agitators, flows in pipes and channels where turbulence is strongly anisotropic and geometry- 
dependent [17-20]. Compared to the few studies of the influence of turbulence in homogeneous and isotropic flow 
conditions (see [7] and [8] among others), the first attempt to assess the effects of flow inhomogeneity and anisotropy 
was recently put forward in the collaborative study of [6]. In this study, direct numerical simulations were performed to 
study the breakage of small inertialess aggregates in different archetypal bounded flows as turbulent channel flow (data 
and analysis provided by this group) and developing boundary layer, comparing results with those of homogeneous 
isotropic turbulence. To emphasize the dependence of breakage on the specific properties of the flow, only brittle 
instantaneous rupture was considered in the simulations. The main findings, relevant to the present discussion, can 
be summarized as follows: (I) Regardless of the flow configuration, the breakage rate decreases when the critical stress 
required to break the aggregate increases; (H) for small values of the critical stress (“weak” aggregates) the breakage 
rate develops a universal power-law scaling that appears to be independent of the flow configuration; (HI) for high 
values of the critical stress (“strong” aggregates) large differences in the breakage rate arise among the different flows 
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and no clear scaling is observed anymore, highlighting the importance of instantaneous and local flow properties. 

Different from [6], in the present paper we examine a more realistic breakage process that results from ductile (non- 
instantaneous) rupture, focusing on the case of turbulent channel flow. We are interested in assessing the influence 
of ductile rupture on the breakage rate, with specific reference to the scaling behavior previously observed for weak 
aggregates [6, 8]. The statistical characterization proposed in this work also provides useful indications about the 
importance of ductile rupture in the context of Large-Eddy Simulations (LES) of industrial particulate systems. 
More specifically, the results discussed here suggest that the specific breakage mechanism has a crucial impact on 
the modelling framework within which LES can be used. As mentioned, shear breakage in anisotropic turbulence is 
controlled by the fluid stresses around the aggregate but these stresses are inevitably filtered out by LES, thus causing 
under-prediction of the breakage rates. This under-prediction is expected to be more evident in the case of ductile 
rupture because the contribution of subgrid fluctuations to the stress history is missing. 

For these purposes, we assume that the breakage process is first activated when the hydrodynamic stress acting on 
the aggregate, referred to as a hereinafter, exceeds a critical value that is characteristic of a given type of aggregate: 
cr > O’er (activation condition, sufficient to produce brittle rupture). In figure 1, which provides a visual rendering of 
the rupture events examined in this study, this condition occurs as soon as the aggregate trespasses the Ocr isosurface 
(point A along the trajectory of the broken aggregate). As long as the condition u > a^r is met the process continues, 
mimicking the situation in which the aggregate is storing energy from the surrounding fluid. The process comes to 
an end when the energy transferred from the fluid to the aggregate, referred to as deformation energy hereinafter and 
defined as: 

E = f e(t\a > CTcrjdt , (1) 

Jo 

with T being the time spent by the aggregate in regions of the flow where a > Ocr and e being the dissipation rate 
of fluid kinetic energy, exceeds the critical breakage value, which is also characteristic of the type of aggregate under 
investigation: E > Ecr (breakage condition). In figure I, this condition occurs at point B inside the acr isosurface. In 
this figure, we also show the trajectory of an aggregate that avoids all regions where a > acr and does not break within 
the time window considered. Note that the acr isosurface is taken at the time of ductile rupture, while aggregate 
trajectories are tracked several time steps backward from this time. To single out the effect of ductile rupture, 
we follow [6] and assume that the stress is cr /r (t/v) ' , where ^ {v) is the dynamic (kinematic) viscosity and 
e = 2vsijSij, with Sij = J {dui/dxj + duj/dxi) the strain rate tensor. Based on these definitions, strong fluctuations 
of e control the fluctuations of the stress and therefore the occurrence of breakage events. In the limit of instantaneous 
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breakage this translates into a picture where an aggregate, once released into the flow, moves through it until the 
local dissipation exceeds a threshold value Ccr that causes its rupture. In this case, the main variable to monitor is 
the time spent by the aggregate in regions of the flow where e < ecr (the so-called exit-time, which will be defined 
formally in the next section). 

We remark here that the proportionality between a and e is strictly valid only for a normal distribution for the local 
shear rate [21, 22], a condition that may be violated in the near-wall region of the channel. In addition, our approach 
still separates the role of turbulence from that of internal stresses, which are neglected. The direct coupling between 
hydrodynamic and internal stresses has been investigated for the first time by [7] for the case of homogeneous and 
isotropic turbulence. Extension of the analysis to anisotropic turbulence is currently under way and will be addressed 
in a subsequent paper. 

II. PHYSICAL PROBLEM AND NUMERICAL METHODOLOGY 

The physical problem considered in this study is the dispersion of tracer aggregates in turbulent channel flow, which 
is the archetypal flow previously analyzed by this group within the benchmark study of [6]. The flow is non-reactive, 
isothermal and incompressible, and the numerical methodology is based on an Eulerian-Lagrangian approach that has 
been used successfully in past investigations of turbulent dispersed flows [23, 24]. The reference geometry consists of 
two infinite flat parallel plates separated by a distance 2h. The origin of the coordinate system is located at the center 
of the channel with the x, y and z axes pointing in the streamwise, spanwise, and wall-normal directions, respectively. 
Periodic boundary conditions are imposed on the fluid velocity field in the homogeneous directions {x and y), no-slip 
boundary conditions are imposed at the walls. The size of the computational domain is LyX Lz = 47rft, x 27r/i x 2h. 
The shear Reynolds number is i?e* = u^^h/v = 150 (Marchioli et ah, 2008), where u* = y'TwJp is the shear velocity 
based on the mean wall shear stress. This value of i?e* was chosen to match that used by [6]. We remark that, based 
on the findings of [25] and those of [26] regarding the statistical distribution of the energy dissipation rate in turbulent 
channel flow up to i?e* = 600, present results are expected to scale up to Reynolds numbers significantly higher than 
i?e* = 150. All variables discussed in this paper are expressed in wall units, obtained using m, and v. 

The flow solver is based on a Fourier-Galerkin pseudo-spectral method that solves for the full Navier-Stokes equa¬ 
tions and thus yields the spatial derivatives required to calculate e along the aggregate trajectory with spectral 
accuracy. Lagrangian tracking is used to calculate the trajectory of each aggregate based on the following equation 
of motion: Xp = UQp, with Xp the aggregate position and U(@p the fluid velocity at Xp. This equation is solved in time 
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using a fourth-order Runge-Kutta scheme, whereas sixth-order Lagrangian polynomials are used to obtain the fluid 
velocity and the fluid velocity derivatives at the instantaneous aggregate position. Further details on the numerical 
methodology can be found in [24, 27]. Breakage measurements are performed by releasing aggregates in two distinct 
regions of the channel: the wall region, which comprises a fluid slab 10 wall unit thick where the viscous stress (rep¬ 
resenting the mean fluid shear) is maximum while the turbulent stress is close to zero; and the center-plane of the 
bulk region, where all wall stress contributions drop to zero and turbulence is closer to homogeneous and isotropic. 
In the following, the two release regions are labeled Vlw and flc, respectively. Within each of these regions, 2 • 10® 
aggregates were released, their trajectories tracked and breakage events detected. 

We remark here that the dissipation of kinetic energy by viscosity in fully-developed turbulence occurs primarily 
at the smallest scales of the flow, namely at scales close of the order of the Kolmogorov scale, tik- In wall-bounded 
turbulence, the mean value of these local dissipation scales can be estimated from the relation ( 2 , t) = /e{z, t)] , 

allowing fluctuations in the scale to be directly connected to variations in e [26]. Due to the non-homogeneity and 
anisotropy of the turbulence [28], dissipation results in local values of e that can be orders of magnitude larger than the 
mean, even for turbulent flows at moderate Reynolds numbers like the one considered in the present study [6]. Such 
high amplitudes are the result of very large velocity gradients that act on the aggregates and ultimately determine 
their breakage. Variations in the smallest scales at which dissipation occurs are reflected in the statistical moments 
of the energy dissipation rate, defined as (e")/{e”), and shown in figure 2 as a function of the distance from the wall, 
z. In this paper, angle brackets denote quantities averaged in time and in the homogeneous directions. Note that, 
because we track many tracer aggregates, the average dissipation seen by aggregates along their Lagrangian trajectory 
is in practice equal to the Eulerian one, which was investigated also by [26] at varying i?e* [29]. Figure 2 shows that 
aggregates are subject to high fluctuations of the kinetic energy dissipation even when they sample the bulk flow 
region. Close to the walls, dissipation attains high mean values while fluctuations, proportional to the Root Mean 
Square (RMS), are intense throughout the channel and correspond to a highly intermittent distribution of e. Figure 2 
also shows that the skewness (n = 3) and flatness (n = 4) factors, S{e) and E(e) respectively, are signihcantly higher 
than both the mean and the RMS, especially in the center of the channel where values differ by roughly three to four 
orders of magnitude. This suggests that breakage events in the bulk of the flow may be caused by (rare) extreme 
energy dissipation events. However, these events are expected to have an effect on instantaneous brittle ruptures more 
than on ductile ruptures, which require the occurrence of events with a certain time persistence. As already observed 
by [26], the higher order moments agree closely with the results in homogeneous isotropic turbulence for much of 
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the channel, exhibiting a universal flow-independent behavior that scales with i?e* and is lost only in the near-wall 
region. This observation can be put in connection with the existence of scaling laws for the breakage rate, observed 
for instance by [6, 8]. Already in the limit of brittle rupture, specific flow properties such as anisotropy and non¬ 
homogeneity have a crucial impact on breakage dynamics since they determine the spatial and temporal distribution 
of fluid stresses (and, therefore, of energy dissipation). In wall-bounded flows such distribution exhibits features 
similar to homogeneous isotropic turbulence in the bulk of the channel, where anisotropy and non-homogeneity are 
not dominant: The behavior of the higher order moments of the energy dissipation shown in figure 1 suggests that 
the breakage process may exhibit universal (or nearly universal) features only in this region. Universality is inevitably 
lost near the wall with important implications for strong aggregates, which can only be broken by extreme dissipation 
events. 


III. RESULTS AND DISCUSSION 

The statistics of the rupture of ductile aggregates are examined by focusing mainly on one observable: the breakage 
rate. Figure 3 shows schematically the procedure we followed to estimate the breakage rate, using the trajectory of 
two sample aggregates labeled as A and B. In figure 3(a), we show the time evolution of the aggregate-to-wall distance 
for A and B, colored using the value (expressed in wall units) of the kinetic energy dissipation along the aggregate 
trajectory: high dissipation is in dark gray (red online), low dissipation is in light gray (blue online). Note that t 
represents the dimensionless time (in wall units) spent by the aggregates within the flow after injection at time to = 0- 
Since we are interested in the effect of flow anisotropy on breakage, the fluid characteristic scale adopted in this study 
is Tf = v/ul- Aggregate A is trapped at the wall and eventually breaks, aggregate B is able to stay in the bulk of the 
flow during the time window of the simulation and escapes breakage. The procedure employed to define the breakage 
rate follows [6, 7] and can be explained with reference to the time history of the dissipation along the trajectory of 
aggregate A, shown in figure 3(b). The aggregate is released at a time to and moves within the flow for a time t = t 
(referred to as exit time hereinafter) after which the local dissipation exceeds for the first time the critical threshold 
ecr, indicated by point A1 in figure 3(b). 

If the aggregate is brittle, then it breaks at point Al: The first crossing of Ccr thus defines the exit time used to 
compute the breakage rate: r = Tbrittie in figure 3(b). The breakage rate for the given threshold is then computed as 
the inverse of the mean exit time, obtained upon ensemble-averaging over many time histories (this averaging being 



7 



Note that, following [6, 7], aggregates are released only at points where the local dissipation is below ecr to avoid the 
occurrence of exit times equal to zero. 

If the aggregate is ductile then it does not break at Al, since e > ecr is just the activation condition. The exit time 
is now computed as the time required to satisfy both the activation condition and the breakage condition, E > Ecr- 
in figure 3(b), the time windows during which both conditions are met encloses all the gray areas, and breakage occurs 
at point A2 yielding an exit time r = TductUe- The breakage rate for the given thresholds of e and E is then computed 


as: 


f (ecr, Ecr) — _ / p , • (3) 

^ductile l,^cr 1 ^cr ) 

Note that, by definition, Ecr depends on ecr- This implies that the statistical distribution of the deformation energy 
made available by the fluid for breakage depends on the specific value of ecr that characterizes the aggregate. Such 
distribution is shown in figure 4, where we plot the Probability Density Function (PDF) of E for aggregates released in 
the channel center, computed according to eq. (1) and considering only aggregates that eventually break. This implies 
that these PDFs provide also evidence of the statistical distribution of breakage events. The curves refer to three 
different values of ecr- £cr = 0.008, corresponding to the case of a weak aggregate for the present flow configuration; 
ecr = 0.7, corresponding to the case of a strong aggregate; and an intermediate case with = 0.12, which can be 
referred to as mild aggregate. As the strength of the aggregate, namely its resistance to breakage, increases, the 
PDF shifts toward higher values of E but exhibits lower peak values (note that PDFs are normalized such that the 
area below each curve is equal to unity). This trend provides a first characterization of the breakage events that are 
typically experienced by the aggregates: very intense but relatively short in time for strong aggregates, less intense 
but more persistent in time for weak aggregates. Note that, since the horizontal axis is plotted in logarithmic scale, 
the PDFs deviate from a Gaussian distribution (as shown clearly in the inset of figure 4, where axes are in linear 
scale): This distribution would be obtained if breakage events were controlled mainly by the early rupture events that 
occur in the vicinity of the aggregate release location (the channel centerline) [6]. Apparently, this is not the case for 
ductile rupture in anisotropic turbulence. 



A. Breakage rates 


To examine further the breakage process, in figures 5 and 6 we show the rates of ductile breakage, obtained according 
to eq. (3) for the two release locations considered in this study: flc and VL\y respectively. These two figures show the 
effect of increasing the critical deformation energy Ecr for different threshold values of the critical energy dissipation 
Ccr and thus extend the findings of [6], which focus on the effect of increasing Ccr (namely the strength of a brittle 
aggregate) in different flow configurations. Results are shown for three different values of the threshold Ecr, sampled 
from the distributions of figure 4 and chosen to span two orders of magnitude: Ecr = 0.04, representing a case of 
low ductility for the present flow configuration; Ecr = 0.4, representing a case of intermediate ductility; Ecr = 4.0, 
representing a case of high ductility. In both figures, breakage rates for brittle aggregates, computed according to 
eq. (2), are also shown for comparison purposes. To ensure convergence of the statistics, the breakage rates reported 
in figures 5 and 6 correspond to a percentage of broken aggregates equal to at least 80% (in the worst case scenario 
of highly ductile rupture, the percentage being above 95% in all other cases). Error bars attached to each profile 
represent the standard deviation from the mean value of the breakage rate, computed using the variance of the exit 
time, = (r^) — (r)^. Error bars are shown to provide an indication of the dispersion of breakage rates around the 
mean value: The higher the dispersion, the lesser the accuracy and predictive capability of correlations based solely 
f {^cr, 

Let us focus first on the results for aggregates released in the center of the channel (figure 5). The breakage rate of 
brittle aggregates (solid curve, taken from [6]) generally decreases with increasing aggregate strength, in agreement 
with the intuitive idea that weak aggregates in wall-bounded flows are broken by turbulent fluctuations faster than 
strong aggregates. Eor small Ccr the breakage rate is known to exhibit a power-law behavior of the type f{ecr) cx 
where x is a flow-dependent scaling exponent: [6] have demonstrated that the value of x for aggregates released in 
the central region of a channel is very similar to that of aggregates released outside a developing boundary layer but 
slightly larger than that of aggregates released near the channel walls or in homogeneous flows. In the case of figure 5, 
the power-law scaling of /(ecr) for brittle aggregates is observed when Ccr < —3 and the best fit is obtained for x — 0.5. 
When ductile aggregates are taken into account (dashed curves), breakage rates change dramatically, especially for 
weak aggregates with low €cr threshold. The values of f {ecr, Ecr) decrease significantly with respect to the case 
of instantaneous breakage, already at low thresholds for the critical deformation energy, Ecr (e.g. Ecr = 0.04). In 
addition, no clear power-law scaling is observed anymore and the breakage rate profiles tend to flatten as the aggregate 
“ductility” increases. As could be expected, the effect of ductile rupture on f {ecr, Ecr) becomes less important for 
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strong aggregates: These must be subject to extremely violent fluid stresses, typical of the intermittent nature of 
small-scale turbulence, to activate the breakage process and thus can store the level of energy required to break 
almost impulsively. As a result, there is just a little increase of the exit time with respect to strong brittle aggregates. 

It is clear from the results of figure 5 that ductile rupture (a process that is of course linked to restructuring phenom¬ 
ena within the aggregate) has a dramatic effect on the frequency with which small aggregates break in wall-bounded 
turbulence. Any predictive model failing to reproduce this feature would inevitably yield strong over-prediction of 
the breakage rates. One example is the exponential model of Kusters [20], which is valid for instantaneous breakage 
only and is based on the dimensional assumption that breakage is ruled by Gaussian kinetic energy dissipation: 

/(Ccr) = 

This classical model predicts a very sharp drop-off at intermediate threshold values of Ccr and a constant breakage 
rate for small threshold values, in strong disagreement to the breakage rate found in our simulations. The discrepancy 
originates both from the simplified assumption of a Gaussian dissipation and from the neglect of ductile rupture. 

The results of figure 5 depend quantitatively on the specihc location chosen to release the aggregates at time to. In 
the center of the channel (release region flc) strong aggregates, no matter if subject to brittle or ductile rupture, are 
mainly broken by the rare extreme excursions of dissipation from the mean, which are caused by intermittency. Most 
of such aggregates must therefore reach the high-dissipation, high-shear regions of the flow near the channel walls to 
undergo breakage. To examine the influence of the release location on breakage rates, in figure 6 we show the behavior 
of /(ccr, Ecr) for aggregates released in the near-wall region Focusing first on the brittle aggregates (solid curve, 

taken from [6]), we observe that the power-law scaling at small values of Scr is followed by a flattening for intermediate 
values of the threshold, which was not observed in figure 5. For the very large threshold values associated to the 
right end of the profile, a drop-off in the breakage rate is observed, representing the case of aggregates that are too 
strong to be broken by the mean shear alone: intense but rare turbulent fluctuations within the near wall region are 
required to overcome the cohesive force of these aggregates [6]. The inclusion of ductile rupture effects (dashed curves) 
produces again a clear decrease of the breakage rates, which vanishes for large values of ecr- Compared to the results 
of figure 5, we observe that the decrease is now almost negligible for aggregates with low ductility (corresponding to 
the Ecr = 0.04 curve in figure 6) and flattening of the profiles is only attained for very high threshold values of the 
deformation energy. We also note that error bars are generally smaller, indicating a lower variability of the statistics: 
this is due to the fact that aggregates are already placed in the high-shear regions of the flow where they preferentially 
break and hence sample a reduced portion of the domain compared to aggregates released in the bulk of the flow. In 
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spite of these quantitative differences, however, the reduction of f{ecr,Ecr) associated with ductile rupture is evident 
independently of the initial position chosen to inject the aggregates into the flow. 

B. Evolution of the number of aggregates 


In the previous section, we pointed out that strong aggregates can move away from the location of release and travel 
towards the high shear regions close to the walls. Obviously, this dynamics has an influence on the breakage process, 
which is also reflected in the time evolution of the number of unbroken aggregates, (t). This quantity can be used 
to derive the following approximation for the breakage rate of brittle aggregates that is valid when breakage is driven 
by homogeneous and temporally uncorrelated stresses [7]: 


/(Ccr) = -- 


dt 


where A^ecr (^) simply linked to the exit time by the relation: 

Ne^rit) 


= 1 -/ PecA'r)dT , 

Jo 


(5) 


( 6 ) 


Nito) Jo 

with N{to) the number of aggregates initially released into the flow and Pe^^ir) the PDF of the exit time for a given 
threshold Ccr- Based on eq. (6), is proportional to the cumulative exit time distribution [7]. 

In figure 7, the evolution of the number of aggregates released in the centerline of the channel is reported. In 
particular the figure shows the behavior of In (t)/N{to)] , considered here because the corresponding slope provides 
a direct estimate of the aggregate breakage rate, as suggested by eq. (5). The different curves refer to a reference 
threshold €cr = 0.008 and to different types of aggregate rupture: brittle (solid line taken from [6]), weakly ductile 
(dashed line), mildly ductile (dotted line) and highly ductile (dash-dotted line). The lowest threshold value for ecr 
was chosen because we know from [6] that in this limit the number of unbroken aggregates decays exponentially as 
~ N(to) exp[—f (ecr) • t], yielding In [A^e^^(t)/A^(to)] — —/(ccr) • t with deviations due only to statistical noise 
(e.g. at late times when the number of aggregates has become very small). The behavior of —/(ccr) • t for the three 
types of ductile rupture examined in figure 7 is represented by the thin solid lines. 

The results of figure 7 show that, in general, the evolution of the number of aggregates follows an exponential decay 
only at short and intermediate times, in good agreement with the linear segments of —/(ccr) • t. Then the decay turns 
into a faster decrease at later times, associated to clear deviations from —f{ecr) • t. These deviations are particularly 
evident for ductile rupture, indicating that simple estimates like the one given in eq. (5) do not provide a reasonable 
approximation of the breakage rate anymore, in agreement with the findings of figure 3. 
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IV. CONCLUSIONS 

In this work, we examine the breakage of small colloidal aggregates in non-homogeneous anisotropic turbulence. 
In particular we focus on the breakage rate of massless aggregates that are subject to ductile rupture caused by the 
hydrodynamic fluid stresses acting on the aggregate. This process is activated when the fluctuating hydrodynamic 
stress generated by the surrounding fluid exceeds a critical value characteristic of a given type of aggregate, cr > acr, 
and ends when the energy given up to the aggregate by the surrounding fluid exceeds the critical breakage value. 
From a physical point of view, the process of ductile breakage comes closer to real applications compared to the case 
of instantaneous rupture considered in previous works [6]. To compute the breakage rate statistics, breakage kinetics 
under a realistic set of assumptions are explored by means of direct numerical simulation of turbulent channel flow 
seeded with a large number of aggregates, modelled as sub-Kolmogorov massless point-particles. Results show that 
the effects associated to ductile rupture are important and lead to strong reductions of the breakage rate with respect 
to instantaneous rupture. The mechanism of ductile breakage thus acts as a low-pass filter for stress-induced events 
that occur at time scales shorter than the characteristic time with which the aggregate responds to deformation. The 
reduction in the breakage rates is evident especially for weak aggregates characterized by small critical stress value and 
no universal scaling can be observed. For strong aggregates characterized by large critical stress value the breakage 
rate is less affected by the specific mechanism leading to rupture because such aggregates can only be disrupted by 
extremely intense stresses and thus store the amount of deformation energy required to break almost impulsively. 

Future investigations will try to evaluate if there is a preferred direction along which breakage takes place. In the 
recent experimental measurements of aggregate breakage in laminar and turbulent shear flows, performed by [30], it 
was found that aggregates tend to break in the direction along which they experience the maximum stretching. Our 
aim is to verify the persistence of this behavior in presence of anisotropic turbulence. 
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FIG. 1: (Color online) Rendering of brittle and ductile rupture in turbulent flow. The trajectory of two different aggregates is 
shown, superimposed onto the isosurface of the critical stress rr = a^r required to produce brittle rupture or activate ductile 
rupture. The broken aggregate trespasses the a^r isosurface at point A (potential brittle rupture) and undergoes ductile rupture 
at point B (where the breakage condition E > Ecr is met). The unbroken aggregate avoids all regions where a > Ocr and does 
not break within the time window considered in this figure. Critical stress isosurface is taken at the time of ductile rupture. 
Aggregate trajectories are tracked several time steps backward from this time. 
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FIG. 2: Statistical moments of the energy dissipation rate seen by the aggregates: mean value, (e); Root Mean Square value, 
RMS{e); Skewness factor, S{e); Flatness factor, F{e). Brackets indicate quantities averaged in time and space over the 
homogeneous flow directions. 
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FIG. 3: (Color online) Time evolution of the aggregate-to-wall distance for two sample aggregates A, subject to breakage, and 
B, not subject to breakage (a); and time evolution of the energy dissipation rate seen by aggregate A along its trajectory. The 
gray areas in panel (b) correspond to time windows during which the aggregate evolves in regions of the flow where the local 
dissipation is above the critical value required to either break the aggregate (in case of brittle aggregate, Tbrittie) or activate 
the breakage process (which, in case of ductile aggregate, is brought to completion at time TducUie)- The trajectories in panel 
(a) are colored based on the instantaneous absolute value of the local energy dissipation rate. Time t is expressed in wall units 
and represents the time spent by the aggregates within the flow after release at time to. 
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FIG. 4: Probability density function of the deformation energy in turbulent channel flow, conditioned to the threshold value of 
ter'- ter = 0.008 (weak aggregate, solid line); ter = 0.12 (mild aggregate, dashed line); ter = 0.7 (strong aggregate, dotted line). 
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FIG. 5: Break-up rate for ductile aggregates released in the center of the channel, f{ecr,Ecr). For comparison purposes, the 
break-up rate of brittle aggregates (solid line, taken from [6]) is shown. Error bars represent the standard deviation from the 
mean value, and were computed using the variance of the exit time, = (r^) — (r)^. 
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FIG. 6: Break-up rate for ductile aggregates released in the near-wall region of the channel, f {ter, Ear)- For comparison 
purposes, the break-up rate of brittle aggregates (solid line, taken from [6]) is shown. Error bars are as in figure 5. 
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FIG. 7: Time evolution of the number of unbroken aggregates (for brittle and ductile rupture of aggregates released in the center 
of the channel). For each type of rupture, the exponential decay predicted by eq. (5) when — N(to) exp[—f (ter) ■ t] 

(thin solid segments) is shown. Solid line taken from [6]. 







